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This  manual  presents  the  complete  ILSS  (revised  ILSLOC)  computer 
program  package.  In  addition  to  including  a thorough  description  of 
the  program  itself  and  a listing  with  comments,  the  manual  contains 
ia  brief  description  of  the  ILS  system  and  antenna  patterns. 
iTo  illustrate  the  program,  a test  case  has  been  created  and  the 
.figures  of  the  case  are/ incorporated  in  the  report.  Program  DYNM 
and  program  ILSPLT  are/included  as  appendixes.  The  ILSPLT,  complete 
with  sample 'graphs , is  a plotting  routine  for  ILSLOC. 

For  a technical  mathematical  analysis  of  the  system,  see  report 
FAA-RD-72-137  (AD754517),  "Instrument  Landing  System  Scattering." 

This  report  revises  in  part  an  earlier  report  FAA-RD-73-76, 
"Users'  Manual  for  ILSLOC:  Simulation  for  Derogation  Effects  on  the 

Localizer  Portion  of  the  Instrument  Landing  System."  The  revisions 
include  the  treatment  of  triangular  scatterers  and  glide  slope 
antenna  systems . . 
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PREFACE 


As  part  of  the  ILS  Perforaance  Prediction  program,  a first 
phase  ILS  Localizer  perforaance  prediction  computer  program  package 
has  been  prepared.  This  package  consists  of  the  computer  program 
and  the  present  document  which  describes  the  capabilities  and 
limitations  of  the  computer  model  as  well  as  the  step  by  step 
running  of  the  computer  program. 

The  computer  program  is  intended  primarily  as  an  aid  in  pre- 
dicting the  perforaance  of  different  ILS  antenna  candidates  for  a 
proposed  runway  instrumentation  or  for  the  upgrading  of  an  already 
instrumented  runway  in  a known  airport  environment.  It  is  also 
intended  to  provide  a relatively  inexpensive  means  by  which  the 
effect  of  proposed  changes  to  an  airport  environment  (addition  of 
terminal  buildings,  hangars,  etc.)  on  ILS  performance  may  be  pre- 
dicted. Another  computer  program  has  been  devised  to  treat  the 
effects .of  terrain  on  glide  slope  performance.* 

This  document  was  prepared  for  the  Transportation  System  Center 

‘W 

(TSC)  by  D.  Newsom  who  was  assigned  as  a full-time  programmer  to 
the  ILS  Performance  Prediction  program.  A.  Watson  and  M.  Scotto 
helped  in  the  writing  and  editing.  The  report  itself  and  the 
attached  computer  program  are  based  on  the  theories  and  analyses 
which  were  developed  by  the  TSC  group  / (Gy/Chin,  L.  Jordan,  D.  Kahn, 
and  S.  Morin).  The  ILS  program  was  spOlisored  by  H.  Butts  of  the 
Systems  Research  and  Development  Service  of  the  Federal  Aviation 
Administration.  ' 

The  present  report  revises  in  part  an  earlier  report,  FAA-RD- 
73-76.  The  revisions  include  the  treatment  of  triangular  scatters 
and  glide  slope  antenna  systems.  The  revised  ILSLOC  program  has 
been  renamed  ILSS-FOR  (Instrument  Landing  System  Scattering- 
Fortran).  The  use  of  the  term  ILSLOC  in  the  body  of  this  report 
refers  to  the  generalized  program,  ILSS. 

4 ^ * 

^S.  Morin  et  al,  ILS  Glide  Slope  Performance  Prediction,  Version  A. 
Report  No.  FAA-RD-74-157  A.  U.S.  Department  of  Transportation, 
Transportation  Systems  Center,  Cambridge  MA  02142,  September 
1974. 
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1.  DEFINITION  OF  INSTRUMENT  LANDING  SYSTEM 


The  ILSLOC  program  has  been  written  to  simulate  certain  air- 
port conditions  which  affect  the  localizer  portion  of  the  instrument 
Landing  System.  The  ILS  is  used  to  provide  signals  for  the  safe 
navigation  of  landing  aircraft  during  periods  of  low  cloud  cover 
and  other  conditions  of  restricted  visual  range.  Separate  systems 
are  used  to  communicate  vertical  and  horizontal  information;  the 
horizontal  system  is  called  the  "localizer." 

This  system  operates  by  the  transmission  of  an  RF  carrier, 
amplitude  modulated  by  two  audio  frequencies,  beamed  to  approaching 
airborne  receivers.  In  an  instrumented  aircraft,  the  localizer 
receiver  serves  to  desK>dulate  the  RP  signal,  amplify  and  isolate 
the  corresponding  audio  signals  and  derive  a signal  to  drive  the 
ILS  horizontal  display  in  the  cockpit.  The  pilot,  by  reading  the 
display,  can  determine  if  he  is  on  course,  to  the  left  of  the 
runway,  or  the  right  of  the  runway.  These  signals  must  be  strong 
enough  to  cover  a radius  of  twenty-five 'miles  around  the  antenna. 

The  directional  information  is  determined  by  the  relative 
strengths  of  the  transmitted  sideband  signals.  The  audio  frequency 
modulations,  which  are  fixed  at  90  Hz  and  ISO  Hz,  are  radiated 
in  different  angular  patterns  with  respect  to  the  runway  centerline 
extended.  The  "course"  is  defined  as  the  locus  of  points  where  the 
amplitudes  of  the  two  modulations  are  equal.  The  display  of  a 
difference  of  the  amplitudes  (90  Hz  and  ISO  Hz)  of  the  sidebands 
is  referred  to  as  the  Course  Deviation  indication.  Thus,  the 
GDI  is  the  pilot's  indication  as  to  what  his  bearing  is  relative 
to  the  center  line  of  the  runway.  The  GDI  is  measured  in  microamps. 
The  actual  course  generated  by  any  particular  ILS  installation  will 
deviate  from  the  ideal  due  to  the  interference  of  spurious  re- 
flections from  buildings  present  in  the  range  of  the  transmitting 
antenna.  The  deviation,  caused  by  these  buildings,  or|  scatterers 
of  the  GDI  from  what  the  receiver  should  read  ideally  at  that 
point  in  space  (e.g.,  on  the  center  of  the  runway  and  GDI  reading 
other  than  0)  is  the  derogation  effect. 
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The  Localizer  systen  transaits  an  asymnetrical  pattern  by 
beaming  a "carrier  plus  sideband"  pattern  and  a "sideband  only" 
pattern,  the  composite  of  wh  .ch  gives  the  desired  effect.  If  a 
specific  localizer  system  ust's  two  antenna  arrays,  four  sets  of 
signals  will  be  transmitted;  if  the  system  uses  a single  antenna 
array,  two  sets  will  be  transmitted. 
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2.  ANTENNA  PATTERNS 


The  proper  angular  variation  of  the  transmitted  90  Hz  and 
the  150  Hz  modulation  is  achieved  by  the  radiation  of  two  independent 
sideband  patterns  by  the  transmitting  antenna  arrays.  Equal 
magnitudes  of  90  Hz  and  150  Hz  modulation  are  transmitted  in  each 
of  these  patterns,  however  with  different  relative  phases.  One 
of  the  patterns  is  symmetrical  with  respect  to  the  prescribed 
course.  An  unmodulated  carrier  wave  is  transmitted  with  the  samd 
pattern  and  the  combination  is  commonly  referred  to  as  the  "car* 
rier  plus  sidebands"  (C  * S)  signal.  The  other  signal  is  trans- 
mitted in  an  "anti -symmetrical"  pattern  and  is  referred  to  as  the 
"sidebands -only"  signal. 

Figure  1 illustrates  how  these  features  are  used  to  obtain 
the  desired  directional  GDI.  The  magnitudes  of  the  C * S and  SO 
sideband  patterns  as  functions  of  angular  deviation  from  the 
course  are  illustrated  in  Figures  la.  The  sideband  amplitude  of 
the  C ^ S pattern  represents  20t  modulation  of  the  carrier  wave 
(or  a "depth  of  modulation"  of  0.2)  at  both  90  Hz  and  150  Hz. 
Considering  the  phases  of  both  modulations  of  the  C S signal  to 
be  positive,  the  relative  phases  and  typical  amplitudes  of  the  two 
SO  modulations  are  as  shown  in  Figures  lb.  The  resultant  90  Hz 
and  150  Hz  modulation  patterns  in  the  total  ILS  signal  are  obtained 
by  algebraically  combining  the  respective  C * S and  SO  sideband 
patterns  (Figures  Ic) . The  evident  consequence  is  that  the  depth 
of  modulation  is  greater  for  90Hz  than  for  150  Hz  to  the  left 
of  the  course  as  seen  from  an  approaching  aircraft,  and  the  op- 
posite is  true  to  the  right  of  the  course.  This  difference  when 
properly  calibrated  in  relation  to  the  total  siodulation  (90  Hz 
* 150  Hz)  reaching  the  aircraft  receiver  give*  the  GDI  as  appears 
in  Figure  Id. 

Since  the  strength  of  C '*'.S  and  SO  signals  fall  off  at  the 
same  rate  with  distance  from  the  transmitting  antenna,  the  GDI 
is  independent  of  range. 
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j ANGLE 


Figure  la  Sideband  Pattern  Magnitude 


Figure  lb  Relative  Amplitudes  and  Phases  In  SO  Pattern 


Figure  Ic  Resultant  Modulation  Patterns 
GDI  I 


Figure  Id  Course  Deviation  Indication  (CDI) 
Figure  1.  Antenna  Patterns  Sketch 
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FAA  standards  for  the  ILS  specify  that  within  a certain  nar- 
row angular  range  about  the  course,  the  CDl  should  be  closely 
proportional  to  the  aircraft's  angular  deviation  from  course.  This 
sector  near  the  ideal  approach  is  termed  the  "course  sector"  and 
usually  extends  between  1-1/2“  and  3"  to  either  side  of  the  runway 
centerline.  The  wider  sectors  on  either  side  of  the  course  sector 
are  called  the  "clearance  sectors."  In  these  sectors,  which  extend 
a minimum  of  35“  from  the  course,  the  CDI  is  required  to  always 
exceed  a certain  minimum  magnitude.  The  presence  of  structures 
in  the  clearance  sectors  which  scatter  spurious  signals  into  the 
course  sector  is  the  primary  cause  of  derogation  of  the  localizer 
CDI.  Such  structures  are  illuminated  by  carrier  and  sideband 
signals.  The  ratios  of  150  Hz  modulation  to  90  Hz  modulation  in 
these  signals  are  determined  by  the  angular  position  of  the; 
structure  with  respect  to  the  runway.  In  general  these  ratios  are 
different  from  those  transmitted  toward  the  aircraft,  due  to  the 
difference  in  angular  position.  The  signals  transmitted  toward 
the  scatterer  will  be  reflected  toward  the  aircraft.  Thus  the 
aircraft  will  receive  the  summations  of  the  direct  and  scattered 
signals.  Since,  in  general,  the  scattered  signals  will  have  im- 
proper ratios  their  effect  is  to  distort  the  CDI.  To  combat  this 
problem  several  new  antenna  systems  have  been  designed.  Two  basic 
systems  are  used:  the  single  antenna,  and  the  "capture  effect 

system." 

The  single  antenna  system  radiates  two  patterns  from  one 
antenna  array.  The  signal  generated  in  the  course  sector  is 
stronger  than  that  generated  in  the  clearance  sector.  However, 
because  of  the  derogation  effects,  the  signals  are  often  not  ac- 
curate enough  to  meet  category  II  or  III  requirements  and  the 
more  accurate  "capture  effect  system"  is  used.  This  system  uses 
one  antenna  array  to  broadcast  a very  narrow,  powerful  beam  in 
the  course  sector.  The  second  antenna  array  broadcasts  a broader 
pattern,  at  a slightly  different  carrier  frequency,  which  covers 
the  clearance  area.  This  system  diminishes  the  derogation  effects 
because  of  the  dual  frequency.  The  term  .'"capture  effect"  has 
been  used  to  describe  this  two-antenna  a/ray  system  because  the 
airplane  receiver  is  "captured"  by  the  stronger  transmission  signal. 
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3.  ILS  SIMULATION  DESCRIPTION 


The  ILS  simulation  program  makes  it  possible  for  airport 

planners  to  determine  what  the  effects  of  potential  airport 

buildings  on  the  ILS  performance  are  going  to  be.  Thus,  for 

example,  if  a new  terminal  or  hotel  is  planned,  the  information 

as  to  size  and  location  of  the  building  can  be  input  to  the  program 

and  the  derogation  effect  of  that  building  can  be  determined. 

Because  the  derogation  effect  of  these  scatterers  is  so  important, 

the  program  can  warn  the  planner  ahead  of  time  to  change  the 

orientation  or  location  of  the  building,  or  it  can  assure  him 

that  the  building  would  not  jeopardize  the  airport's  current 

FAA  rating.  » 

* ! 

The  output  of  this  program  is  a magnetic  tape  of  values  of 
the  GDI.  Graphs  are  generated  by  a plotting  routine  (ising  the  ^ 
values  derived  from  the  ILSLOC  program)  to  show  the  CDI  in  micro- 
amperes, along  a flight  path,  for  the  scattering  surfaces  input. 
These  generated  graphs  would  serve  the  same  purpose  as  the  FAA 
strip  charts  which  are  generated  for  a certifying  flight.  The 
simulation  graph  differs  from  the  actual  recorded  measurements 
due  to  limitations  of  the  program  which  will  be  explained  later 
in  the  text. 

The  ILSLOC  program  simulates:  transmission  from  the  various 

types  of  localizer  antenna  systems;  the  trajectory  of  an  aircraft 
flight  over  which  the  CDI  is  to  be  determined;  and  the  scattering 
from  rectangular  and  cylindrical  surfaces.  The  program  permits 
various  simulated  flight  paths. 

The  program  is  not  an  exact  simulation  of  the  certifying 
flight,  due  to  certain  simplifying  assumptions  which  were  made. 
These  assumptions  include: 

a.  A flat  perfectly  conducting  ground  plane 

b.  Perfectly  conducting  reflectors 
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c.  |Far-field  scattering--  all  scattering  fron  a surface 

is  assumed  independent  of  all  other  surfaces;  thus, 
multiple  reflections  from  walls  and  near-field 
interactions  are  ignored 

d.  A noise-free  environment 

e.  Relative  field  strengths--  the  absolute  field  strengths 
involved  are  not  calculated;  thus  while  we  can  calculate 
the  GDI's  in  microamperes,  we  do  not  ascertain  the 
absolute  electric-field  intensities,  and 

f.  An  idealized  ILS  receiver  model. 

In  addition  to  these  assumptions  the  approximations  of  the 
scatterer  can  lose  accuracy  when  the  dimensions  approach  less 
than  a few  wavelengths.  Since  the  program  determines  the  scat- 
tering from  a surface  independently  from  all  other  scatterers,  the 
shadowing  of  one  structure  on  another  is  not  included.  Thus  if 
one  building  is  between  the  antenna  system  and  another  building, 
it  will  shield  the  second  one  from  some  of  all  of  the  ILS  signal. 
The  amount  of  energy  reaching  the  second  building  will  depend  upon 
diffraction  effects  which  are,  in  general,  too  complicated  to 
analyze.  It  may  be  noted,  however,  that  diffraction  effects 
themselves  are  included  as  part  of  the  physical  optics  approxima- 
tion used.*  By  using  rule  of  thumb  approximations  the 
analyst  can  determine  roughly  how  much  power  will  reach  the 
second  building.  If  the  level  is  small  the  building  may  be 
ignored  completely.  If  on  the  other  hand  the  power  level  is 
large  then  the  structure  should  probably  be  included  as  though 
there  was  no  shielding  effect.  This  will  give  a conservative 
GDI  estimate  (i.e.,  larger  derogation  than  actual),  but  this 
will  serve  for  most  purposes.  If  the  situation  is  critical, 
that  is  near  category  limits,  then  other  means  of  analysis  must 
be  used.  ' 

Chin,  G.,  L.  Jordon,  D.  Kahn,  and  S.  Morin, 

TSC,  "Instrument  Landing  System  Scattering," 

FAA-RD-72-137,  AD754517  (Dec.  1972). 
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4,  TEST  CASE  FOR  THE  ILSLOC  COMPUTER  PROGRAM 


To  illustrate  how  the  computer  program  is  operated  a very 
simple  test  case  (with  only  2 scatters)  has  been  created  and  run. 

For  this  simulated  airport  the  program  computed  the  course  width 
as  4.01  degrees.  Both  antenna  arrays  were  set  at  an  elevation 
of  13  feet  above  the  ground  plane.  The  clearance  antenna  array 
was  used  as  the  origin  for  the  coordinate  system.  An  80  x 100  x 60 
ft  hangar  and  75  x 110  ft  cylinder  were  placed  on  opposite  sides 
of  the  9,350  ft  runway.  In  this  case  the  threshold  is  10,000  ft  ' ^ 
from  the  course  antenna.  (See  illustration--Figure  2.)  Based  on 
the  size  and  locations  of  these  two  buildings,  the  model  pre- 
dicted the  CDI  on  the  runway  centerline  and  for  a clearance  run  at 
10,000  ft  range. 

Using,  this  model  for  input  values,  the  following  section 
presents  a detailed  follow  through  of  the  main  program  steps. 

The  Mode  Card 

The  first  input  is  the  mode  card.  This  card  .contains  informa- 
tion on  the  type  of  localizer  antenna  used,  the  frequency  of  the 
ILS,  the  length  of  the  runway,  and  the  height  of  the  antenna.  The 
mode  card  format  is  shown  immediately  following  Figure  2.  In  order 
to  use  the  mode  card,  it  is  important  that  the  user  understand  the 
coordinate  system  used.  The  x-axis  is  along  the  centerline  of  the 
runway,  the  threshold  being  in  the  positive  direction.  The  z-axis 
is  vertical,  positive  z being  in  the  up  direction.  The  y-axis 
completes  a right-hand  coordinate  system;  so  that  when  one  is 
standing  at  the  origin  facing  in  the  x-direction,  positive  y is 
to  the  left.  The  origin  is  used  as  a reference  to  define  the 
location  of  scatterers,  antenna  system  components,  and  flight 
path  sample  points.  The  antennas  are  located  along  the  x-axis, 
they  need  not  be  at  the  origin.  As  in  our  test  case,  it  is 
usually  convenient  to  place  the  course  antenna  at  the  origin. 
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Model  Card  Format 


Col.  Symbol 

1-2  Mode 


3-4  lET 

11-20  FKQ. 

21-30  XTH 

31-40  ZA(1) 

41-50  ZA{2) 

51-60  SLP 


Use 

■ 1 
- 2 

- 3 

- 4 

- 5 

- 6 

- 7 

- 8 

--1 

--2 

--3 

--4 

-13 

-14 

-15 

- 0 

- 1 


(V-Ring) 

(8-Loop) 

(Wave  Guide) 

Not  Used 

(Measured  Pattern) 
(Measured  Capture 
Effect  Patterns)  i 

(Theoretical  Patterns)  ^ 
(Theoretical  Capture  f 
Effect  Patterns) 

(V-Ring  Clearance 
(8-Loop  Clearance) 
(Waveguide  Clearance)  \ 
(Measured  Clearance 
Patterns) 

Null  Reference 
Sideband  Reference  ] 

N-Array  (Capture  Effect)* 


J 


Indicates 

Localizer 

Antenna 

Type 


Glide  Slope 
Antennas 


(Half-wave  dipole)  Antenna 

(30“  width  pattern  as  from  \ Element 
double  wave  reflector  used] 
with  glideslope  antennas  J Pattern 


Frequency  of  ILS  in  MHz 

Distance  from  the  origin  to 
the  threshold  of  the  runway, 
in  feet.  This  number  is  used  ' <> 

for  both  flight  path  orienta- 
tion and  for  course  width  de- 
termination. The  distance  is 
given  in  feet. 


There  is  always  a non-zero 
antenna  height,  and  xt  is 
input  here. 


This  will  be  the  clearance 
antenna  height  if  a two- 
antenna  system  is  used. 


Slope  of  runway  in  degrees 
used  to  adjust  the  flight 
path  at  threshold  for 
runway  slope. 
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Modes  1,  2,  and  3 provide  for  standard  localizer  antenna 
types.  These  antennas  are  predetermined,  the  only  variable  being 
course  width,  the  adjustment  of  which  is  controlled  by  the  course 
width  card. 

When  any  antenna  type  other  than  mode  1,  2,  or  3 is  used, 
additional  antenna  description  card.^  must  be  included.  Mode  5 
permits  the  input  of  a measured  pattern  for  special  cases  on 
theoretical  studies.  When  this  mode  is  selected  additional  pattern 
cards  are  required.  One  pattern  card  must  be  used  for  each  measure 
ment.  The  angles  must  be  given  in  ascending  order.  A maximum  of 

50  measurements  may  be  given;  if  less  than  50  cards  are  used  a 
termination  card  with  an  angle  greater  than  360  degrees  must  be 
inserted. 

Format  of  Pattern  Card(s) 


Col.  Symbol  Use 

1-10  ANG  Angle  of  measurement,  in  degrees 


11-20  AFPP  Amplitude  of  sideband  only  pattern, 

in  relative  units 

21-30  AGPP  Amplitude  of  carrier  plus  sideband  pattern, 

in  relative  units. 


Mode  7 allows  the  generation  of  a theoretical  array  pattern 
from  assumed  element  contributions.  The  antenna  is  to  be  a linear 
array  of  elements  with  identical  radiation  patterns.  Each  element 
has  an  arbitrary  magnitude  and  phase  for  both  carrier  plus  side- 
band and  sideband  only  currents.  The  arrays  are  assumed  to  be 
alligned  parallel  to  the  y-axis.  All  elements  have  the  same  height, 
as  given  in  the  mode  card.  All  elements  have  the  same  x-coordinate 
as  given  on  the  course  width  card.  The  y-coordinate,  in  wavelengths, 
is  given  for  each  element  on  the  element  description  card.  There 
must  be  one  card  for  each  element  in  the  array,  to  a maximum  of  26 
elements.  The  element  description  card  has  the  following  format: 


11 


Col. 

Symbol 

Oser 

1-10 

DT 

Element  displacement  in  the  y-direction 
given  in  wavelengths 

11-20 

CT 

Carrier  plus  sideband  amplitude,  in 
relative  units 

21-30 

PC 

Carrier  plus  sideband  phase,  in 
degrees 

31-40 

ST 

Sideband  only  relative  amplitude 

41-50 

PS 

Sideband  only  phase,  in  degrees. 

The  phase  of  the  sideband  only  currents  is  ideally  in  quad- 
rature to  the  carrier  plus  the  sideband  currents.  This  90-degree 
shift  is  added  by  the  program.  Thus  a "PS"  inputted  as  zero  degrees 
is  internally  converted  to  90  degrees  out  of  phase  with  the  sideband 
portion  of  the  carrier  plus  sideband.  To  indicate  termination  when 
there  are  less  than  26  elements  used,  an  element  card  is  placed  with 
a carrier  plus  sideband  phase  value  (PC)  of  more  than  SOO. 

The  next  step  for  this  mode  must  be  the  input  of  the  horizontal 
radiation  pattern  for  the  individual  element.  This  pattern  will  be 
used  for  each  of  the  elements  previously  described.  The  input  is 
the  relative  signal  strength  measured  every  10“  starting  at  0 and 
proceeding  until  180“.  This  is  total  of  19  amplitudes;  the  values 
are  read  in,  in  records  of  8F10.4  format,  for  a total  of  3 records. 
This  gives  the  pattern  for  angles  from  0 to  180“  and  since  the 
pattern  is  assumed  to  be  symmetric  the  value  for  the  negative  angle 
will  be  the  same  as  a positive  one  of  equal  magnitude. 

There  are  two  methods  of  inputing  capture  effect  system 
descriptions.  The  most  general  way  is  to  input  each  antenna 
separately.  When  using  this  method  the  clearance  antenna  must  be 
input  first.  This  input  will  follow  the  same  steps  as  a single 
antenna  system  except  that  the  mode  number  will  be  a negative. 

The  negative  mode  card  and  the  pattern  or  element  cards  (if  any) 
must  be  followed  by  another  mode  card.  This  mode  for  the  course 
antenna  must  be  positive,  and  followed  by  the  necessary  pattern  or 
element  cards. 


There  are  two  cases  for  the  second  Method  of  inputing  antenna 
descriptions.  The  first  case  is  used  if  both  course  and  clearance 
antennas  are  to  be  given  as  measured  patterns;  a single  mode  6 
card  is  used  followed  by  two  sets  of  pattern  cards:  the  first  set 
is  for  the  course  antenna  and  the  second  set  for  the  clearance 
antenna.  In  the  second  case,  for  a capture  effect  system  which 
uses  two  theoretical  array  antennas,  a mode  8 is  used.  This  card 
is  followed  by  the  course  antenna  element  description  cards  and 
the  element  radiation  cards;  a second  set  of  array  description 
cards  is  used  in  the  clearance  antenna. 

In  the  above  localizer  antenna  cases,  lET  has  no  effect  on 
the  simulated  individual  element  patterns  and  may  be  input  as 
zero. 

FRQ  is  set  to  the  frequency  (in  MHz)  of  the  carrier  for  the 
antennas  system. 

XTH  is  the  distance  (in  feet)  from  the  orgin  to  the  runway 
threshold.  The  flight  path  is  set  to  cross  the  threshold  at  an 
altitude  of  ZUP  (as  given  on  flight  path  card). 

ZA(1)  (course)  and  ZA(2)  (clearance)  are  the  heights  in  feet 
of  the  antennae. 

SLP  is  the  slope  of  the  runway  in  degrees.  It  is  used  with  XTH 
in  setting  up  the  flight  path.  The  ground  plane  assumed  for  the 
signal  scattering  is  not  tilted. 

Modes  13,  14,  and  15  are  used  for  glide  slope  antennas.  Al- 
though this  program  is  intended  for  localizer  simulation,  it  may 
be  also  used  to  study  the  effects  of  builuings  on  the  glide  slope 
system.  The  simulation  will  assume  a perfect  flat  horizontal  and 
infinite  ground  plane.  If  a glide  slope  antenna  is  chosen  on  the 
mode  card,  the  next  card  must  be  as  follows: 


Col 

Symbol 

Use 

1-10 

YA 

Antenna  Offset  (in  feet) 

11-20 

TGS 

Glide  path  angle  (in  degrees) 

13 


YA  is  the  antennae  offset  (Y-coordinate)  in  feet.  Positive  is  to 
the  left  from  the  origin  facing  the  threshold.  If  the  magnitude 
of  YA  is  less  than  300,  YA  will  be  defaulted  to  1500,  the  sign 
depending  on  the  sign  of  YA  input.  TGS  is  the  intended  glide  path 
angle  in  degrees. 

The  measured  pattern  of  a capture  effect  localizer  is  used 


in  our  test  case: 
Mode  Card: 

Col.  1-2 

6 

11-20 

110. 

21-30 

10000. 

31-40 

13. 

41-50 

13. 

Pattern  Cards:  See  attached  Figure  3 for  test  case  listing. 

The  antenna  description  cards  are  followed  by  the  course 


width  card. 

The  format  for  this 

card  is : 

Col. 

Symbol 

Use 

1-10 

XXA(l) 

Course  array  x-coordinate , 
in  feet 

i 

11-20 

XXA(2) 

Clearance  array  x-coordinate, 
in  feet  , 

31-40 

CW 

Course  width  in  degrees 

41-50 

CLS 

Clearance  signal  strength 

relative  to  the  course  signal. 


If  CW  is  greater  than  3”  this  value  is  used  as  the  course 
width  and  the  signal  strengths  of  the  course  antenna  are  auto- 
matically adjusted  to  produce  this  value. 

If  CW  is  less  than  3*  the  course  width  will  be  set  to  the 
FAA  specification  for  a threshold  to  antenna  distance,  given  by 
XTH,  and  the  signal  levels  will  be  set  accordingly. 

CLS  is  the  ratio  of  clearance  signal  strength  to  course 
signal  strength. 
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-45. 

-.012 

0.006 

-42. 

-.020 

0.014 

-40. 

-.014 

0.020 

-38. 

0.000 

0.020 

-35. 

0.O1B 

O.ooo 

-32. 

0.008 

-.025 

-30. 

-.010 

-.020 

-28. 

-.011 

0.000 

• 

-27. 

-.008 

0.010 

-26. 

0.000 

0.017 

-25. 

0.011 

0.019 

-23. 

O.020 

0.000 

-20. 

0.000 

-.039 

-19. 

-.010 

-.041 

-18. 

-.015 

-.035 

-16. 

0.000 

0.000 

-14. 

0.016 

0.024 

-13. 

0.015 

0.035 

-12. 

0.000 

0.050 

-9. 

-.180 

0.140 

-5. 

-.535 

0.535 

-4. 

-.535 

0.660 

-1. 

-.165 

0.996 

* I. 

n. 

0.000 

1.000 

1. 

0.165 

0.996 

4. 

0.535 

O.f.60 

- 

5. 

0.535 

0.535 

9. 

0.180 

0.14O 

12. 

0.000 

0.050 

1 

13. 

-.015 

0.035 

1 

14. 

-.016 

0.024 

16. 

0.000 

O.OOO 

1 

18. 

0.015 

-.035 

19. 

0.010 

-.<'41 

20. 

0.000 

-.039 

1 

23. 

-.020 

0.0<'o 

1 

25. 

-.011 

0.010 

k 

26. 

0.000 

0.017 

i 

27. 

0,008 

o.on 

28. 

0.011 

o.or-''  1 

i 

30. 

0.010 

-.020 

i 

i « 

32. 

-.008 

-.025 

35. 

-.018 

o.ooc  ' 

38. 

o.ooo 

O.O20 

40. 

0.014 

o.r'2'' 

42. 

0.O20 

0.01 

43. 

0.O12 

O.Ol  p 

1000. 

1 

Fifur*  3. 

Pattern 

Card  Test  Casa  Listing 

IS 

-60, 

O.rtOO 

0,000 

-55. 

-.085 

0,018 

-54. 

-.096 

0.019 

-51. 

-.145 

0,008 

-50, 

-.160 

0.002 

-49, 

-.175 

0,005 

-45, 

-.245 

0.050 

-33. 

-.411 

0.400 

-32. 

-.414 

0.430 

-30. 

-.426 

0.475 

-27. 

-.464 

0.497 

-26. 

-.475 

0.499 

-25. 

-.440 

0.497 

-22. 

-.545 

0.486 

-21. 

-.565 

0.485 

-20, 

-.585 

0.486 

-14. 

-.602 

0.490 

-15. 

-.676 

0.540 

-14. 

-.680 

0t560 

-13. 

-.680 

0.585 

-12. 

-.675 

0.620 

-9. 

-.610 

0.730 

-2. 

-.160 

0.980 

0. 

0,000 

l.OOO 

2. 

0,160 

0.980 

4. 

0.610 

0.730 

12. 

0,675 

0.620 

13. 

0,680 

0.585 

14. 

0.680 

0.560 

15. 

0.676 

0.540 

19. 

0.602 

0.490 

20. 

0,585 

0.486 

21. 

0,565 

0.485 

22. 

0.545 

0.486 

25. 

0,490 

0.497 

26. 

0,475 

0.499 

27. 

0,464 

0.497 

^0. 

0,426 

0,475 

32. 

0.414 

0.430 

33. 

0.411 

0.400 

45. 

0.245 

0.050 

49. 

0,175 

0.005 

50, 

0.160 

0.002 

51. 

0.145 

0.008 

54. 

0.096 

0.019 

55. 

0.085 

0.018 

60, 

0,000 

0.000 

1000, 

Figure  3.  Pattern  Card  Teat  Caae  Listing  (Cont'd) 
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The  test  case  course  width  card  would  read: 


1-10 

0 

11-20 

-200 

31-40 

0 0 

41-50 

0.315 

The  label  card  follows  the  course  width  card.  This  card  is 
put  on  the  output  tape  ahead  of  the  CDI  records  for  this  flight. 
It  serves  as  an  identifying  record  and  is  the  label  placed  on  the 
graph.  Columns  180  are  used.  In  our  test  case  this  card  reads: 
THIS  IS  A DEMONSTRATION  CASE  OF  STRAIGHT  LINE  FLIGHT. 

The  program  calculates  the  COI  at  a point  in  space:  for 

convenience,  the  program  will  permit  calculation  for  a series 
of  points.  This  set  of  points  represents  samples  of  a simulated 
flight  path. 

The  program  allows  two  types  of  flight  paths.  A straight 
line  flight  and  a circular  orbit.  The  flight  path  card  has  one 
of  the  following  formats: 


Straight  Line  Flight 


Col. 

Symbol 

Use 

1-10 

XMIN 

Starting  distance  from  origin, 
in  feet 

11-20 

XMAX 

Ending  distance  from  origin,  in 
feet 

21-30 

DXR 

Spacing  between  sample  points, 
in  feet 

31-40 

PHIR 

Angle  of  approach,  in  degrees 

41-50 

PSIR 

Glide  angle,  in  degrees 

61-70 

ZUP 

Height  of  aircraft  at  threshold. 

in  feet. 


XMIN  is  the  x*coordinate  of  the  starting  location  of  the 
aircraft  and  XMAX  is  the  x-coordinate  of  the  ending  location. 
The  sample  points  are  spaced  along  a straight  line  so  that  the 
difference  in  x*coordinates  between  successive  samples  is  DXR. 
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The  sign  of  the  DXR  will  be  set  by  the  program  so  that  the 
flight  goes  from  XMIN  to  XMAX  regardless  of  flight  direction. 

If  the  DXR  value  would  require  more  than  500  points  the  program 
will  adjust  the  magnitude  of  DXR  to  give  only  500  points.  In 
some  cases  a flight  will  require  more  than  500  points.  If  this 
is  necessary  the  flight  must  be  broken  up  into  smaller  segments 
of  not  more  than  500  points  each.  The  procedure  for  doing  this 
is  explained  in  the  control  card  section.  The  flight  path  is 
oriented  in  space  so  that  an  extension  of  the  path  crosses  the 
threshold  at  the  altitude  of  ZUP  and  intersects  the  z-axis.  PHIR 
is  the  angle  between  the  flight  path  and  the  vertical  plane  through 
the  runway  centerline.  It  is  zero  for  a flight  path  along  the 
centerline  of  the  runway  and  is  positive  for  an  incoming  flight 
(XMIN  greater  than  XMAX)  with  decreasing  y-displacement.  PSIR 
is  the  glide  angle  between  the  flight  path  and  the  horizontal 
plane.  It  is  zero  for  level  flight  and  positive  for  a normal 
landing  approach.  The  flight  path  is  a straight  line  as  de- 
scribed above  except  when  the  x-component  is  less  than  XTH,  that 
is  if  the  aircraft  is  on  the  antenna  side  of  the  threshold.  In 
that  case  the  aircraft  altitude  will  be  set  up  to  ZUP. 

Thus  the  values  used  in  the  test  case  would  read; 

40000 
20000 
-40 
0 

2 5 
50 

The  arc  flight  is  a series  of  points  at  a constant  height 
of  ZUP  and  at  a constant  horizontal  distance  from  origin  of  R. 

MIND  is  the  starting  angle  for  the  arc,  that  is,  the  line  of 
sight  from  the  origin  to  the  point  makes  a horizontal  angle  of 
MIND  degree  with  the  x-axis.  The  sample  points  are  spaced  at 
equal  angles  of  DXR  until  the  termination  angle  of  MIND  is 
reached.  As  in  the  straight  line  flight  the  sign  of  DXR  will 
be  adjusted  appropriately.  Likewise  the  magnitude  of  DXR  will 
be  set  to  yield  not  more  than  500  points.  Column  74  must  be 
set  to  1 to  indicate  a circular  arc. 


Col.  1-10 
11-20 
21-30 
31-40 
41-50 
51-60 
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r 


Circular  Orbit  Case 


: 


Col. 

Symbol 

Use 

1-10 

MIND 

Starting  angle,  in  degrees 

11-20 

MAXD 

Ending  angle,  in  degrees 

21-30 

DXR 

Angular  spacing  between  samples, 
in  degrees 

51-60 

R 

Radius  of  orbit,  in  feet 

61-70 

ZUP 

Height  of  orbit,  in  feet 

74 

ICF 

Must  be  set  to  1 to  indicate 
orbit  case. 

Following  the  flight 
in  the  following  format: 

path  card  must  be  the  velocity  card 

Col. 

Symbol 

Use 

1-10 

VEL 

Velocity  of  aircraft,  in  feet/sec. 
This  is  used  for  the  Doppler  Effect 
on  the  receiver.  The  sign  of  the 

velocity  will  be  made  to  agree  with 
the  directional  motion  from  DXR.  Test 
case  assumes  velocity  of  200  ft/sec. 

At  this  point  we  have  described  the  antenna  system  and  the 
trajectory  of  the  aircraft;  the  derogating  surfaces  in  proximity 
to  the  ILS  must  now  be  described.  The  program  will  simulate 
scattering  from  rectangular  or  cylindrical  surfaces.  We  will  now 
describe  the  method  of  inputting  scatterers  to  simulate  derogating 
structures . 

The  next  card  describes  either  the  scatterer(s)  or  output  and 
control.  The  usage  is  determined  by  the  value  of  the  ID  field 
in  columns  1 to  2.  An  ID  of  -1 , 1*  or  2 is  used  for  scatterers, 
while  the  other  values  are  used  for  control. 
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Col. 

Symbol 

Use 

1-2 

ID 

Must  be  1 for  rectangle 

3-8 

XW(1) 

x-coordinate  of  reference  point, 
in  feet 

« 

9-14 

XW(2) 

y-coordinate 

15-20 

XW(3) 

z-coordinate 

26-30 

ALPHA 

Angle  between  base  and  x-axis,  in 
degrees 

31-35 

DELTA 

Angle  of  tilt,  in  degrees 

36-45 

WW 

Width  of  rectangle,  in  feet 

46-55 

HW 

Height  along  rectangle,  in  feet. 

The  scatterer  is  a rectangle  with  the  reference  point  at  the 
middle  of  the  base.  The  rectangle  is  assumed  to  be  of  infinite 
conductivity  and  zero  thickness.  It  also  has  only  one  side.  This 
can  be  thought  of  as  the  front  surface  of  a metal  wall.  A wall 
with  zero  x- , y-,  and  z^coordinates  and  an  alpha  of  zero  is  located 
at  the  origin  with  surface  of  the  wall  facing  in  the  negative  y 
direction  (Figure  4,  case  I).  A positive  increase  in  alpha  rotates 
the  wall  about  the  z-axis  in  a counterclockwise  direction  when 
viewed  from  above.  Thus  an  alpha  of  90  degrees  faces  the  wall 
in  the  positive  x-direction  (Figure  4,  case  II).  Alpha  is  the 
angle  between  the  vertical  projection  of  the  base  of  the  wall  in 
the  xy-plane  and  the  x-axis,  measured  in  degrees.  Delta  is  the 
angle  between  the  surface  of  the  wall  and  the  vertical  direction, 
in  degrees.  A delta  of  zero  is  a wall  perpendicular  to  the  ground 
and  a decrease  in  delta  rotates  the  wall  about  the  baseline  in  a 
direction  so  that  a delta  of  minus  ninety  is  a horizontal  wall 
facing  down  (Figure  4,  case  III).  WW  is  the  width,  in  feet,  of 
the  wall  measured  along  its  base  and  HW  is  the  height  measured 
along  the  surface  at  right  angles  to  the  base.  If  the  wall  is 
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oriented  in  such  a fashion  that  the  line  of  sight  from  the  antenna 
to  the  wall  passes  through  the  back  and  not  the  front  of  the  wall, 
the  program  will  ignore  the  wall  in  the  simulation. 

An  ID  of  -1  is  used  with  the  above  format  to  describe  a 
negative  wall.  This  ID  is  used,  for  example,  to  create  a wall 
with  a rectangular  hole  in  it.  The  entire  surface  is  used;  the 
hole  is  then  subtracted  by  inputing  a second  card  with  an  ID  of 
-1  and  the  size,  location,  and  orientation  of  the  hole. 

An  ID  of  2 is  used  for  a cylindrical  scatterer  with  the 


following  format: 

Col. 

Symbol 

Use 

1-2 

ID 

Must 

be 

a 2 

3-8 

XW(1) 

9-14 

XW(2) 

n 

y 1 
1 

coordinates  of  the 
reference  point,  in  feet 

15-20 

XW(3) 

z-  , 

i 

- 

36-45 

WW 

Diameter 

of  cylinder,  in  feet 

46-55 

HW 

1 Height  of  cylinder,  in  feet. 

The  reference  point  is  located  at  the  base  of  the  cylinder  on 
the  axis  of  rotation  of  the  cylinder.  The  diameter  is  WW  feet,  with 
the  base  parallel  to  the  xy-plane  at  an  altitude  of  XW(3)  feet.  The 
cylinder  extends  upward  for  HW  feet  with  the  axis  of  rotation  in  the 
vertical  direction.  The  cylinder  is  assumed  to  have  infinite 
conductivity. 


t 
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An  ID  of  3 or  -3  is  used  for  triangular  scatters  with  the 
following  format: 


Col. 

Symbol 

Use 

1-2 

ID 

Must  be 

3 or  -3 

3-8 

xw(l) 

X- 

9-14 

xw(2) 

Y- 

Coordinates  of  the  reference 
point,  in  feet. 

15-20 

xw(3) 

z- 

26-30 

Alpha 

Angle 

between  base  and  x-axis,  in  degrees 

31-35 

Delta 

Angle 

of  tilt,  in  degrees 

36-45 

WW 

Width  of  base  of  triangle,  in  feet 

46-55 

HW 

Height 

along  side  of  triangle,  in  feet. 

The 

variables  have  the 

same  use  as  for  the  rectangular  scat- 

terer,  with  the  exception  of  HW  4 WW.  The  magnitudes  of  WN  S HW 
determine  the  size  of  the  triangles,  the  signs  of  HW  4 WW  are 
used  to  determine  the  orientation  of  the  hypotenuse.  The 
convention  is  as  follows: 


Sign  of  Sign  of 

Triangle  Orientation  HW  WW 

1^ 

- - 
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If  the  size  of  the  triangle  exceeds  the  limits  imposed  by  the 
Fresnel  approximation  the  scatterer  will  be  omitted  and  an  error 
message  given  in  the  output.  If  this  happens  and  one  wishes  to 
include  scattering  from  this  surface,  the  triangle  must  be  broken 
up  into  triangular  and  rectangular  pieces,  for  example: 


The  values  for  IH  and  IV  will  indicate  the  number  of  pieces  horizon- 
tally and  vertically  the  triangle  must  be  broken  up  into. 
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After  an  ID  of  1,-1,  2,  3,  or  -3,  the  program  will  calculate 
the  electric  field  at  the  surface  of  the  scatterer.  This  will  be 
calculated  from  the  signal  from  the  transmission  antenna  array 
and  from  the  ground  reflection  of  the  transmitted  signal.  Then, 
for  each  receiver  point  along  the  flight  path,  the  program  will 
calculate  the  electric  field  at  that  location  from  the  scattered 
signal:  from  both  the  scatterer  and  reflected  from  the  ground. 

Thus,  the  signal  is  received  from  four  paths:  (a)  transmission 

antenna  to  scatterer  to  receiver,  (b)  antenna  to  ground  to  scat- 
terer to  receiver,  (c)  antenna  to  scatterer  to  ground  to  receiver, 
and  (d)  antenna  to  ground  to  scatterer  to  ground  to  receiver.  This 
signal  is  decomposed  into  complex  components  induced  in  the 
receiving  antenna  at  the  different  carrier  and  sideband  fre- 
quencies. The  program  then  loops  back  to  read  in  another  ID 
card,  permitting  the  summation  of  the  effects  of  many  scatters. 

This  allows  the  simulation  of  complex  structures  by  breaking 
them  up  into  cylinders  and  rectangles. 

In  the  test  case,  we  have  only  inputed  three  scattering 
surfaces.  This  was  done  because  only  two  sides  of  the  hangar 
and  the  cylinder  are  illuminated.  The  values  for  the  scatterer 
cards  read: 

Col.  First  card  Second  card  Third  card 

1-21  1 2 


3-8 

6000 

5950 

7500 

9-14 

1100 

1130 

-1000 

15-20 

0 

0 

0 

26-30 

10 

-80 

0 

31-35 

36-45 

100 

60 

75 

46-55 

80 

80 

110 

25 


After  all  the  scatters  have  been  input,  a control  card  is 
inserted  to  terminate  the  run.  The  control  card  format  is: 

Col . Symbol  Use 

1-2  ID  not  -1,  1,  or  2. 

When  a control  card  is  read  in,  the  program  will  add  the  direct, 
and  ground  reflected,  signal  from  the  transmission  antenna  to  the 
scattered  signal  summations,  thus  giving  the  total  received  signal. 
The  program  then  calculates  the  CDI  that  would  be  seen  at  each  re- 
ceiver point,  and  outputs  the  label,  a header  record  describing 
the  flight  path  and  the  values  of  the  CDI  on  output  tape.  If  the 
ID  is  equal  to  zero  the  program  also  outputs  additional  records  for 
the  strengths  of  sideband  and  carrier  signals  from  course  and 
clearance  (if  any)  antenna  arrays.  The  field  summations  are  then 
cleared  for  the  next  run. 

The  program,  having  finished  the  previous  run,  now  proceeds 
with  the  next  input.  The  next  run  is  generated  by  looping  back 
to  a point  in  the  input  stream,  determined  by  the  value  on  the 
control  card. 

Once  an  input  sequence  has  begun  the  inputs  following  in  the 
standard  order  must  be  given.  The  user  must  also  keep  in  mind 
that  all  values  on  cards  given  before  that  entry  point,  in  the 
previous  run  are  still  in  effect.  The  following  order  is 
standard : 

NODE  CARD 

(measured  pattern  for  modes  5 and  6 or  current 
description  for  modes  7 and  8) 

(second  mode  card  and  patterns  of  currents  if 
first  mode  was  negative) 

COURSE  WIDTH  CARD 

LABEL  CARD 

FLIGHT  PATH  CARD 

VELOCITY  CARD 

(set  of  scatterer  cards) 

CONTROL  CARD. 
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The  value  of  the  ID  on  the  control  card  guides  the  looping 
in  the  following  manner: 

Value  of  ID  Next  card  to  be  read  in 

0 MODE 


3-10 

SCATTERER 

11-15 

LABEL 

16-20 

MODE 

21-50 

COURSE  WIDTH 

>50 

WILL  CAUSE  THE  PROGRAM  TO 
TERMINATE  AFTER  OUTPUTTING 
THE  LAST  CDI. 

The  looping  permits  the  repetition  of  a run  with  changes  in 
some  or  all  of  the  variables.  For  example,  ID  values  3 through  10 
permit  a run  with  the  same  antenna  system  and  flight  path  as  the 
previous  case,  but  with  a new  set  of  scatterer  inputs. 

ID  values  11  to  15  permit  a new  flight  path  description  and 
scatterer  set  to  be  input.  This  looping  method  can  also  be  used 
for  flights  that  would  require  more  than  500  points.  For  reliable 
simulation,  the  spacing  between  receiver  points  (DXR)  should  be 
small  enough  so  that  the  change  in  CDI  between  successive  points 
is  not  more  than  ~20t  of  the  peak  value.  Thus  for  long  flights  the 
flight  path  must  be  broken  up  into  shorter  segments.  If  the  number 
of  segments  of  this  path  does  not  exceed  4,  the  plotting  program  will 
connect  them  on  a single  graph.  The  control  for  this  joining  is  the 
ID  number.  If  the  flight  path  finishes  with  an  ID  of  11  to  13,  the 
graph  of  the  next  flight  will  continue  the  line  of  the  graph.  A 
long  flight  may  be  broken  up  into  as  many  as  four  segments:  with 
three  segments  terminating  in  11  to  13  and  a fourth,  and  final  seg- 
ment, terminating  in  14  or  15.  The  flight  segments  must  appear  in 
the  order  in  which  they  are  to  be  flown,  so  that  the  XMIN  of  one 
section  is  the  XMAX  of  the  previous  section.  For  each  segment 
the  programmer  must  re -input  the  same  scatterers.  If  only  one 
segment  is  to  be  plotted  the  control  card  should  read  14  or  15. 
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id's  16  through  20  start  inputing  at  the  mode  card,  thus 
allowing  a completely  new  run. 

An  ID  of  21  through  SO  uses  the  same  antenna  description,  but 
starts  the  inputing  at  the  course  width  card.  This  permits  the 
course  width,  clearance  strength  and  antenna  location  to  be 
varied. 

The  program  is  terminated  after  an  ID  greater  than  50  is  en- 
countered. The  direct  signal  will  be  added,  and  the  CDI  will  be 
outputed  before  the  program  stops.  The  program  will  also  stop 
if  an  end  of  file  is  encountered  while  the  program  is  attempting 
to  read  any  input  card,  or  if  certain  of  the  variables  are  of  im- 
proper value.  In  these  cases  the  program  terminates  immediately, 
without  outputing  the  last  case. 

The  input  of  the  test  case  flight  path  was  done  in  four 
segments.  The  first  segment  is  from  40,000  to  20,000  feet,  the 
second  segment  is  from  20,000  to  12,500  ft,  the  third  segment  is 
from  12,500  to  11,000  ft  and  the  last  is  from  11,000  to  10,000  ft. 
An  additional  case  for  a simulated  clearance  flight  by  a circular 
orbit  has  also  been  included.  The  input  cards  for  these  test  case 
flights  are  shown  in  Figure  5. 
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THIS 

IS  A DEMONSTRATION  CASE 

OF  straight 

LINE  FLIGHT 

AOOftO, 

20000, 

-AO. 

2.5 

200, 

16000. 

1100. 

10. 

100. 

80. 

15950. 

1130. 

-80. 

60. 

80. 

27500. 

-1000,  0, 

0. 

0.  75. 

110. 

13 

THIS 

IS  A DEMONSTRATION  CASE 

OF  STRAIGHT 

LINE  FLIGHT 

20000. 

12500, 

-15. 

2.5 

200, 

16000, 

1100. 

10. 

100. 

80. 

15950, 

1 130. 

-80, 

60. 

80. 

275on, 

-1000,  0, 

0. 

0.  75, 

110. 

13 

THIS 

IS  A DEMONSTRATION  CASE 

OF  STRAIGHT 

LINE  FLIGHT 

12500. 

11000, 

-3. 

2.5 

290, 

16900, 

1100. 

10. 

100. 

80. 

$0. 


30. 


50. 


15950.  1130. 

-80, 

27500.  -1000,  0, 

0. 

0 

13 

THIS  IS  A DEMONSTRATION 

CASE 

OF 

11000.  10000, 

2. 

290, 

16990.  1100, 

10, 

15950.  1130, 

-80. 

27500.  -loon,  0, 

0, 

0 

15 

60. 

80. 

75. 

110. 

STRAIGHT 

LINE  FLIGHT 

2.5 

50. 

100. 

80. 

60. 

80. 

75, 

110. 

THI 

1«0. 

200, 

16000. 

15950. 

27500. 


S IS  ORBIT  CASE  WITH  SIGNAL 
IBO,  0,72 


1100. 

11?0, 

-1000.  0, 


10. 

-80. 

0.  0. 


STENGTHS 


100. 

60. 

75. 


10000. 


80. 

80. 

110. 


50. 


Figure  S.  Flight  Case  Inputs 
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The  ILSLOC  prograa  calculates  the  CDI  at  each  point  in 
space;  this  CDI  includes  the  Doppler  effects  from  the  velocity 
of  the  aircraft.  In  the  siaulation,  the  receiver  system  is 
assuned  to  generate  the  CDI  value  instantaneously.  In  the  real 
case,  the  inertia  of  the  electrical  and  mechanical  portions  of 
the  system  limit  the  rate  of  change  of  the  CDI.  Thus  the  real 
observed  CDI  appears  to  have  been  low-pass  filtered  from  the 
instantaneous  CDI. 

The  program  DYNM  takes  the  output  tape  generated  by  program 
ILSLOC  and  converts  it  to  observed  CDI  by  simulating  the  effect 
of  a low-pass  filter.  The  variable  TAD  is  the  time  constant  of 
the  effective  filter.* 


Note: When  a flight  path  has  been  segmented,  the  low-pass 
filter  will  operate  continuously  over  the  entire 
flight  path. 


C THIS  PROGRAM  SIMULATES  THE  EFFECT  OF  THE  MECHANICAL  AND  ELECTRICAL 
C INERTIA  OF  THE  ILS  RECEIVER  ON  THE  CDI.  THIS  EFFECT  IS  EQUIVALENT 
C TO  A simple  R-C  low  PASS  FILTER*  THE  VARIABLE  TAU  IS  THE  TIME 
C constant  of  the  effective  filter,  a typical  value  is  .4  SECONDS* 

C THE  INPUT  tape  IS  ON  UNIT  11.  THE  OUTPUT  ON  UNIT  12* 


DIMENSION  XV ( 10 ) ,OEF< 501 ) .MEMO! 14) 

LOGICAL  FOF 
DATA  ILBL/4HDYNM/ 
data  TAU/0,4/ 

IF(EOF(lin  GO  TO  4 

1 IT*?' 

DELC-0. 

2 READ(ll.lOOO)  MEMO.XY.ID.NC.ICF 
WRITE(6,1003)  MEMO.XY.ID.NC.ICF 
DEFIC*ABS(XY(9)/XY(5)/TAU) 

IR*IFIX(XY(10)+.l) 

READ! 11.1001)  (DEF(I).I>1.IR) 

IFiIT  *E0*  0)  CEF2-DEFI1) 

IT*1 

DO  3 1*1. IR 
CEF2-CEF2+OELC 
DELC*(DEF( I)-CEF2)*DEFK 

3 DEF(I)-CEF2 

MEM0(13)*ILBL  . 

WRITFI 12.1000)  MEMO.XY.ID.NC.ICF 
WR1TEI12.1001)  (DEF( I) .I*1.IR) 

IF(ID  *GT*  13)  GO  TO  1 
IF(  ID  *EQ*  0)  GO  TO  1 
GO  TO  2 

4 REWIND  11 
END  FILF  12 
REWIND  12 
CALL  EXIT 

1000  FORMAT(13A6.A2./.1X.7F18.9./.3F18*9.I10.10X.2I10) 

1001  FORMAT(7E15.8) 

1003  F0RMAT(1X.13A6.A?./.1X.7F18*9./.3F18.9.I10.10X.2I10) 
STOP 
END 
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APPENDIX  C 

ILSPLT  PLOTTING  ROUTINE 

4 


1 


l(  This  program  has  been  written  to  generate  graphs  of  the 

static  and  dynamic  GDI's.  It  was  written  on  the  IBM  7094  using 
the  CALCOMP  plotting  subroutines. 

The  first  input  card  has  the  following  format: 


Gol. 

Symbol 

Use 

1-2 

NL 

Number  of  lines  per  graph 

3-4 

NGRFS 

Number  of  graphs 

5-7 

NTAPE  (1) 

Input  logical  unit  no.  for 
line 

first 

8-10 

NTAPE  (2) 

Input  logical  unit  no.  for 
line 

second 

11-13 

NTAPE  (3) 

Input  logical  unit  no.  for 

third 

line. 

NL  permits  the  overlaying  of  two  or  more  GDI  or  signal 
strength  graphs  for  comparison  purposes.  The  scaling  will  be 
set  by  the  first  graph,  and  the  successive  overlays  will  be 
plotted  to  the  same  scale.  A maximum  of  three  lines  per  graph 
will  be  allowed. 

NGRFS  sets  the  maximum  number  of  graphs  to  be  drawn.  Each 
graph  will  have  the  same  number  of  overlays. 

NTAPE  (i)  gives  the  logical  unit  number  used  for  the  input 
of  the  ith  line  on  each  graph.  If  the  value  of  NTAPE  is  negative 
then  its  absolute  value  will  be  used  as  its  logical  unit  number 
and  the  tape  will  be  rewound  before  input. 

The  second  input  card  defines  the  scaling  used  for  the  graph 
(or  graphs)  described  above.  It  has  the  following  format: 
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Col. 

Symbol 

Use 

1-10 

XSC 

Horizontal  scale  in  ft/in.  or  deg/in. 

11-20 

DELX 

Tick  mark  spacing  in  ft.  or  deg. 

21-30 

YMAX 

Maximum  y-value  on  vertical  scale 

31-40 

YMIN 

Minimum  y-value  on  vertical  scale 

41-50 

DELY 

Tick  mark  spacing  on  vertical  spacing 
in  microamps  for  CDI  or  relative  units. 

The  horizontal  axis  is  drawn  in 'either  feet  or  degrees  per 
inch  as  specified  by  XSC.  The  ticl:  mark  spacing  along  the  axis 
is  determined  by  OELX.  The  length  of  the  axis  will  be  adjusted 
to  the  shortest  length  with  an  integral  number  of  tick  marks 
that  will  cover  the  domain  required  by  the  input  data.  When  a 
flight  path  has  been  segmented  it  is  treated  as  a single  line 
on  the  graph. 

YMAX,  YMIN  define  the  range  of  the  plotted  variable:  CDI 

or  relative  signal  strength.  The  y-axis  has  a fixed  length  of 
seven  inches.  If  DELY  does  not  integrally  divide  the  range,  DELY 
will  be  adjusted  to  yield  an  integer.  When  the  range  (YMAX-YNIN) 
is  zero,  the  program  will  automatically  scale  the  range  to  the 
largest  scale  that  will  include  the  data  in  the  length  of  the 
axis . 

When  multiple  graphs  are  plotted,  each  graph  is  scaled  in- 
dependently. 

After  all  NGRFS  graphs  have  been  drawn,  the  program  will  loop 
back  to  the  beginning  and  attempt  to  read  in  a new  NL  card.  This 
allows  many  graphs  to  be  drawn.  If  the  user  wishes  to  replot 
data  using  different  scales  or  overlaid  with  different  sets  of 
data,  he  may  use  the  negative  NTAPE  to  rewind  the  input  tape. 

The  program  will  terminate  after  reaching  an  end-of-file 
on  the  card  input  unit. 

The  vertical  scale  on  the  graph  is  always  labeled  "micro- 
amperes." This  is  valid  only  for  CDI  graphs.  All  others  are 
in  relative  units  and  this  labeling  should  be  ignored. 
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COMHCN/TEST/XMINtDXfttNrOT»NP 
LOGICAL  EOF 
DIMENSION  IBUF(IOOO) 

DMENSICN  NTAPEOl  tNEP40(L4ltM(14l 
EQUIVALENCE  (Ml 1| .NEMOl 111 

CCMMCN  /PDF/  OFIZOOOIfXLENtNSTEPSf iDEFtlOENTtOXf 101 iNPTSIlOl 

COMMON  /PAINT/  ML i XSC tOELX *VMAX« VMIN. OELY t ICF 

CALL  PLOTS! IBUF, 10001 

CALL  FL0TI0.0t-12.t-3l 

CALL  FACTOR  (0.41 

ILBL«1 

60  CONTINUE 

IF(E0F(5li  GO  TO  55 
RFAD(5tl00l  NLtNGRFStNTAPE 
WRITE(6tl00l  NLtNGRFStNTAPE 
IFINGRFS.LE.CI  NGRFS-3 
100  F3RMAT(212t3T3l 

00  401  l«ltNL 

IF(NTAPF(II.GE.0I  CO  TO  401 

NTAPEdl—NTAPEdI 

NU«NTAPEd) 

RCWINO  NU 
401  CONTINUE 

REA0(5tlCl I XSCtOFLXtVMAXtVMINtOELY 
WRITE(6t 1911  XSCtOELXtYNAXfYMINtOELY 
101  FORMAT(BFIO.O) 

TEMP«AMINl(VPINtYMAXI 
YMAX-AMAXl I YMlNt  YMAXI 
YMIN-TEMP 
TEMP«VPAX-YM1N 

1F(TEMP  .NE.  0.1  0ELY-TEMP/(FL0ATdFIX(TEMP/0ELY».5l  II 
NPLT  - I 
NP  ■ 1 

1 * 1 
N1  - 1 
NTOT  ■ 0 

10  NU  - NTAPE(NPI 

IF(EOF(NUII  GO  TO  50 

REA0(NUtS001  M.XOtOXRtXYtlOtlOEFtlOENTtICF 
IFdCF  .NE.  01  ICF-1 

WRITE (6t6COI  MEMOtXOtOXRtXYtlOtlOEFt lOENTtlCF 

1F( ILBL  .NE.  II  GO  TO  70 

ILBL«0 

CALL  SYMBOL! 0.t0.t.l4tMEM0t90.«80l 
CALL  PLOT!3.t0.t-3l 
70  CONTINUE 

IR  >IFIX!  XT*’.!! 

NTOT  - NTOT  ♦ IR 
IFd.EQ.ll  XMIN  • XO 
900  FORMAT! 13A6tA2t/t/t3F18.9f4I 101 
600  F0RMAT!2Xt 13A6,A2t/t 3F 18.9*41 101 
501  FCRMATdE19.8l 

902  FORMAT! lXt7E19.8l 

REAO!NUt901l!0FIJIt J«NltNTOTI 
WRITE!6t902l  !0F!JI t J«NltNT0TI 


99 


WRlTE(6t 1000)  XMIN«IR»Nl»NTOTtNPtl 
1000  F0RMAT|F10.3»5110) 

NPTSm  > IR 
OXf  I } • OXR 

1F(  10  «GT.  13  ) GO  TO  40 
If  no  .EO.  0)  GO  TO  40 
N1  * N1  ♦ IR 
I ■ I ♦ 1 
GO  TO  10 
11  NL  « NP 
40  CO.^TUUE 
NSTEPS  ■ I 

IFINP.GT.l)  GO  TO  41 
CALL  GRAPH2(0i 
GO  TO  42 

41  CALL  GRAPH2(1) 

42  CONTINUE 
N1  « 1 
1 » 1 
NTOT  ■ 0 

IF(NP.EO.NL)  GO  TO  45 
NP  ■ NP  ♦ 1 
GO  TO  10 
45  NP  - 1 

CALL  PL0TIXLEN»7.f-12.i-3l 

NPLT  ■ NPLT  ♦ 1 

ILSL*! 

IFINPLT.GT.NGRFS  ) GO  TO  60 
GO  TO  10 
50  CONTINUE 

IFINTOT.GT.O)  GO  TO  11 
CALL  PLOT  <XLEN*7.t-12.#-3l 
GO  TO  60 
55  CONTINUE 

CALL  PLOTIO.tO. »999) 

00  400  I«ltNL 
NU«NTAPE<1) 

4C0  RE  NINO  NU 
STOP 
ENO 


SUBROUTINE  CRAP H2 1 ITU 
OINENSICN  XLABU) 

COMMON /TE ST /XO* DEL TAXt NOE LTAtNP 

DATA  XLAB/24H0ISTANCEtFT.  DEGREES  / 

OIMENSICN  TVPEIBI 
DIMENSION  XI 3). NCI  31 

COMMON  /PDF/  OFI2000ltXLENtNSTEPStIOEFtlOENTtOX(10l,NPTSI10l 
COMMON  /PRINT/  NL »XSC »DELX »YMAX» VMIN» OELYt ICF 
DATA  X />S..S.t5./ 

DATA  NC  /l*5tA/ 

IFIITL  .NE.  C)  CO  TO  1 
ELX-DELX 

I^IOELTAX.LT.O.)  ELX  - -ABSfOELXI 
RANGE-0. 

DO  11  I-ltNSTEPS 

11  RANGE«RANGE»FLOAT INPTS II 11*0X1 1) 

TIX-1FIXIRANGE/ELX«^.9I 

7 XLFN  « ABSIELX/XSC*TIX) 

IFIXLEN  .ST.  40.)  GO  TO  9 
IFIXLEN  .GT.  5.)  GO  TO  6 

9 XSC-A8SIRANGE/20.) 

XLEN«ABSIELX/XSC*TIX) 

WR1TFI6»8)  XSC 

8 F0RMATI25H  AXIS  OUT  OF  RANGE  SCALE- «E 12. 5. 8H  FT. /IN.  /) 

6 CONTINUE 

XMAX-TIX*ELX»XO 
XMIN  > APINlIXOtXMAX) 

XMAX  - AMAXllXOtXMAX) 

NO  - 2 
PMR  « 0. 

CALL  PLCTIO. tl.5«-3) 

AMIN-YMIN 

AMAX-YPAX 

IFIYMAX  .EO.  YMIN)  CALL  SCLAXI 7. t OF ,NOELTA » AMAX* AMI N, OELY* NOt PWR 1 
CALL  AXIS3IO.tO.»AMAX»AMIN»OELV«7.»12HMICROAMPEREStl2»NOtPMRtOELN) 
YSC  - OELN 
IXLAB«2*ICF*1 
IXSC— 1 

IFIABSICLX)  .LT.  10.)  IXSC-1 

CALL  AXlS3IO.t0.tXMAXtXMlNfELX»-XLEN*XLABIlXLAB) *12  >IXSC*0. 
.fOELN) 

XSC  - OELN 

XT  • XLEN/2.  - 2. 

IFIAM1N*AMAX.GT.0.)  GO  TO  2 
1F|  AMIN  .EQ.  0.)  CO  TO  2 
2ER0-I0.-AMIN/10.**PNR)/YSC 
CALL  PLOT  I0.»2ER0»3) 

CALL  PL0TIXLEN*ZER0»2) 

2 CONTINUE 
1 CONTINUE 

XI-0. 

IFIOELTAX  .LT.  0.)  XI-XMAX-XMIN 
J-1 

03  9 I-1*NSTEPS 
OELTAX  > 0X1 T) 
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NX-NPTSf n 

IFII  ,LT.  NSTEPSI  NX-NX*1 
VM-APIN/10.**PWR 

CALL  XCLINECXItOELTAXtOPIJ) iNXtO. tXSCtVM,YSC«NCINP) ) 
J«J*NPTSm 

XI-Xl«-OXm*FLOATINPTStI)> 

5 COMTINUE 
RETURN 
END 
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SUBROUTINE  XCLINEIXI tOX *Y »N»XM,DEUtVNt OELY»NC ) 
DIMENSION  Yl  U»IPEN<4) 

REAL  L(4,4)»LL(4) 

DATA  IPEN/2f3t2»3/ 

OAT  A L/«3f*lt*3i«lt*St  3P»  OS* • 3 *3P* ltalt*OS*«lt*  05/ 
X > XI 

2 IC  ■ NC  - I 

XPl  ■ (X-XM)/DELX 
YPl-IYm-VN)/OELY 
CALL  PLCT(XP1*YP1«3) 

IF(IC.LE.O)  GO  TO  1000 
IF(IC.GT.4)  IC  « 4 
K-1 
“ 1-2 

X - X ♦ OX 
XP2  « (X-XM)/DELX 
YP2»(Y<2»-YM)/0ELY 
1 LL(K)«L(K,IC) 

10  0IFFX-XP2-XP1 
DIPFY-YP2-YPI 

DIS»S0RT(0IFFX*0IFFX4DIFFY*0IFFYI 
IF(DIS.GT.LL(K))60  TO  100 
CALL  PLCT(XP2,YP2,1PEN(K|) 

XP1-XP2 

YP1-YP2 

1-I*1 

IF( I.GT.NIRETURN 

X - X ♦ OX 

XP2  « <X-XM>/OELX 

YP2-IYII »-YN»/DELY 

LL(K)-LL(K)-01S 

GO  TO  10 

IOC  RATIO-DIS/LL(K) 

XP1-XP1»0IFFX/RATI0 
YP1-YP140IFFY/RAT10 
CALL  PL0T(XP1*YP1*1PEN(K|) 

K-K  + 1 

1FCK.EQ.S)K>1 
GO  TO  1 

1000  00  SO  I>2*N 
X - X ♦ OX 
XPl  - (X-XMI/OELX 
YPl- IY(1 )-YHI/OELY 
SC  CALL  PLCTIXP1«YP1«2I 
RETURN 
END 
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SUBROUTINE  SCLAXUINCH* VAR«NtVMAXtVHINtOELTA»NOtEXP) 
DIMENSION  VARll) 

AXLEN  - AINCH 
VMAX  - VARIll 
VMIN  ■ VARin 
DO  40  I«2»N 

VNAX  « AMAXKVHAX.VARm) 

40  VMIN  « AMINKVMINtVARlIM 
NO  « 0 

NE  « 0 
M - 2 

TOTAL  « VMAX  - VMIN 

DETERMINE  EXPONENT  AND  INCREMENT/ INCH 
VM  * AHAXl(ABSIVMAX)tABSIVNIN)) 

Iff VMAX*VMIN}  6*5t7 
7 VAV  « ABSf  VMAX1-VNINI/2. 

DELTA  ■ TOTAL/AXLEN 

If (TOTAL. 6T.0.. AND. TOTAL/VM.LT. .751  GO  TO  4 
IF(VMAX.EQ.VM)  VMIN-0. 

IFCVMIN.FO.-VH)  VMAX-0. 

GO  TO  5 

6 AXLEN  » AXLEN*VM/TOTAL 
5 DELTA  - VM/AXLEN 
VAV  ■ VM/2. 

TEST  FOR  VAV  BETWEEN  .01  AND  1000. 

4 IF(VAV.LE.1.E-11I  GO  TO  21 
IFIVAV  - .01)  3>10»1 

41  IFIVAV  • 1.)  BflO.lO 

1 IFIVAV  - 1000.)  10.2t2 

VAV  GE  lOCC. 

2 IFINE.EO.O)  VAV  - VM 
VAV  « VAV/1000. 

NE  ■ NE  - 3 
GO  TO  1 

VAV  LT  1. 

3 VAV  - VAV«1000. 

NE  - NE  «■  3 

GO  TO  41 

oetermine  decimal  places  in  delta 

10  IFI0ELTA.LT.VM/1.E4)  GO  TO  21 
DELTA  - 0ELTA410.44NE 

11  IFIOELTA  • 1.)  12fl9tl3 

12  DELTA  « OELTAMIO. 

NO  > NO  ♦ 1 

GO  TO  11 

13  IFIOELTA  • 10.)  15«8il4 

14  DELTA  > OELTA/10. 

NO  ■ NO  • 1 

GO  TO  13 

DELTA  NOW  BETWEEN  1 AND  10 

19  IFIOELTA  • 9.1  1A*1T«17 

16  IFIOELTA  - 2.)  19«lBtlB 

17  DELTA  ■ 9./1C.64IN04NE) 

GO  TO  20 
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18  DELTA  - 2./10.**IN0«NEI 
H • 5 

GO  TO  20 
8 NO  - NO  - 1 

19  DELTA  - l./10.**IN0»NE) 

C KESET  VNIN  lEIRSTV)  fOU  AXIS 

20  AK  ■ VNIN/OELTA  ♦ .01 

K « ( IFIXt AK t/NI*N 

IFIVM1N.lt. 0.)  K-K-M 

VNIN  « OELTA*FLOAT«KI 

NOIV  - IVMAX  - VMINI/OELTA  ♦ .9 

IF(FL0AT(N01V).GT.A1NCH*2.)  DELTA*0ELTA*ANAXl(2.»FL0ATIN)/2.l 
1F(N0.LE.0I  NO  - -1 

21  EXP  » NF 

W9ITE(6»1002)  VMAXtVMINtOELTAtNOtNE 
AETUAN 

1002  P09NAT<lHf»,3E13.3»3I7// ) 

END 


lOS 


SUaROUTINC  AXIS3IX0»Y0»VMAX»VMIN*DELVtRINCHtSC0tNCRtN0ECtPWRtVSC) 

FACTOR  ■ lO.WKR 

AMIN  - VN1N«FACT0R 

AMAX  • VMAX«FACT0R 

OELX  - ABS(OELV)*FACTOR 

oiMENSicN  acorn 

HT  ■ .15 

W1«0. 

M2«0. 

M3  - 0. 

NEXP  - 0 
NCH-IABS(NCRI 
SFtPMR.NE.O.I  NEXP  > 6 
CINCH«ABS(AINCH) 

1F(| VNAX-VMIN|/AMAX1IVMAX*-VMIN).LT.1.E-6I  go  to  50 
IFIfAMAX-AMIN)/(DFLX>l.E>81.GT.3.*CINCH)  OELX  > ( AMAX-AMINI /Cl NCH 
IFCOEU.GT.AMAX-AMINI  OELX  - AMAX  - AMIN 
IFINCR.LT.O)  M3  « 1. 

NUM-(AMAX-AMIN)/DELX«1.9 

ANC>CINCH/FL0AT(NUM-1) 

IFIAINCH.lt. C.IGO  TO  5 

W2«l. 

GO  TO  10 
5 Mi*l. 

19  CALL  FL0T(X0»Y0.3) 

VSC  « OELX/FACTOR/ANC 
ANUM«AM1N>0ELX 

X-0. 

Y«0. 

XM>0. 

OFF  - .05 
on  40  1<1,NUP 
ANJM>ANUH>OELX 
II«0 

25  IFIABSIANUM)/10.**n.LT. 1.160  TO  20 
1I«1I*1 
GO  TO  25 

20  IF(ANUP.LT.C.II1«I1»1 
iFuesf anumi.lt. 1.)  II«lt«>l 
IM0RE*N0EC«>1 
I1«II«IM0RE 

IF<IF|XIN1|41.EQ.1I  HT  • AMINllHT  »ANC/FLOATf  1 U2II 
HL  ■ AMAX1(.12*1.2*HTI 
CENTER  • FLOAT! 11 l•HT/ll•♦Ml) 

XC  ■ X - CENTER  - M2*. 15 

IFIXC.LT.XMl  XM  ■ XC 

1F(M2*H3.GT.0.1  XC  - .15 

IFIABSlXCl.GT.ABSIXMl)  XM  ■ XC 

YC  • Y - M1*«MT  ♦ .15  - N3*(HT*.31I  - M2*0EF 

CALL  PLOTlXOAXtYO^Yt21 

CALL  PLOTIXO«XP.l*N2tYO*Y«.l*Ml*3) 

CALL  PLOT I XO»X-. 1* M2  »YO*V-. l*Nl • 2 1 
CALL  NUMBER (XOPXCtY0^YCfHT»ANUf«»0.*N0ECl 
CALL  PLCTIX0*X»Y0Py»3l 
X*X^ANC*lil 
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V>V^ANC*W2 
40  CONTIMUE 

SST  « ICINCH  - FLOAT(NCH«-NEXP)*HL)/2. 
mWi.EO.l.)  XN  *>  -XM 

XXC  - W1*(X0  ♦ BSTl  ♦ W2*IX0  ♦ XP  - OFF  «-  M3*t2.*0FF>HL II 
YVC  ■ ♦ YC  - 1.5PHL  ♦ W3*<HT  ♦ 2.*HLll  ♦ W2*«Y0^SSTI 

CALL  SYMBrLIXXCtYYCtHL«BC0t90.*W2tNCHI 
IFIPMR.EO.O.I  RETUBN 

CALL  SVMBCL(999.t999..HL»5H  * 10t9C.*W2tSI 
X - 999.  ♦ IXXC-.669HL-999.19W2 
V - 999.  ♦ IYYCP.66*ML-999.l*Bl 
CALL  NUMBER ( X, Yt. 79*HL tPWR *90. *M2«>1I 
RETURN 

5C  VSC  « I VMAX-VMIN*l.E-6/FACT0RI/CllilCH 
WRITE(6*1U00I 

1000  FnRMATUH0*27HINSUFFtCIENT  RANGE  FOR  AXIS  I 
RETURN 
END 
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SIMULATED  TEST  FLIGHT  SHOWING  EFFECTS 


PATTERN  - SIDEBAND  ONLY  SHOWING  SCATTERERS 


